The production of K^K pairs in proton-proton collisions at 2.83 GeV 
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Differential and total cross sections for the pp — >■ ppK^ K~ reaction have been measured at a 
proton beam energy of 2.83 GeV using the COSY-ANKE magnetic spectrometer. Detailed model 
descriptions fitted to a variety of one-dimensional distributions permit the separation of the pp — ^ 
pp(j> cross section from that of non-0 production. The differential spectra show that higher partial 
waves represent the majority of the pp — >■ ppcf) total cross section at an excess energy of 76 MeV, 
whose energy dependence would then seem to require some s-wave (fyp enhancement near threshold. 
The non-(f) data can be described in terms of the combined effects of two-body final state interactions 
using the same effective scattering parameters determined from lower energy data. 

PACS numbers: 13.75.-n, 14.40.Be, 25.40.Ep, 13.75.Jz 



I. INTRODUCTION 

The phenomenological description of strangeness pro- 
duction in nucleon-nucleon collisions near threshold is 
complicated for a variety of reasons and these add to 
its general interest. The large mass changes involved 
are necessarily associated with short-range phenomena 
and therefore stress the importance of heavy meson ex- 
change. Whether such exchanges are mainly of strange 
or non-strange nature is still an open question. 

It is known that the scattering lengths in the Ap and 
K^p systems are both very large and these will distort 
any spectra. Furthermore, important transitions, such 
as ^ Ap, mean that several of the final channels are 
strongly coupled. It is also possible that, due to such 
a coupling, the production of a heavy hyperon in pp — > 
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i4'+pA(1405) might influence kaon pair production in the 
pp — > ppK~^K~ reaction fH. The only hope of being 
able to disentangle such effects is through having detailed 
experimental spectra in different kinematic variables. 

In addition to explicit strangeness production one has 
to consider also hidden strangeness, such as that resid- 
ing in the meson which, in the quark model, is mainly 
composed of ss pairs. The pp — pp<j) cross section might 
therefore be influenced by other strangeness production 
channels that are important at this energy. We know 
from the Okubo-Zweig-Iizuka (OZI) rule that the (p 
production rate should be much lower than that of the 
Lj meson in the pp — >■ ppuj channel. Some of the widely 
observed violations of this rule, such as for example in 
the 4>/uj ratio measured in pp annihilation Q, might be 
understood if there were significant strangeness compo- 
nents in the nucleon. There are, however, alternative 
explanations in terms of modified meson exchange mod- 
els [1, The production of the meson in proton- 
nucleus collisions will clearly depend on the more ele- 
mentary pp — >■ pp4) cross section and so this will also be 
an important ingredient in the understanding of proton- 
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induced nuclear transparency measurements [6|. 

The production of the (j) meson in pp collisions has been 
studied in several theoretical papers within a simple one- 
meson-exchange model 3, including both mesonic and 
nucleonic currents components IsMlll. a nd also contri- 
butions from nucleon resonances [12l - ll4| . However, the 
rather limited published data set is not sufficient to pro- 
vide strong constraints on the different models. It is the 
aim of the present paper to present detailed measure- 
ments of the pp — ?> ppK'^ K~ reaction at a beam energy 
of 2.83 GeV, where the cross sections for the production 
of the 4> meson is cleanly separated from the non-0. 

We have previously published data on the pp — >■ 
ppK^ K~ reaction [15, 16] at excess energies with respect 
to the threshold, — ^~{2mp+m^)c^ , of 18.5, 34.5, 
and 76 MeV, where \/s is the total center-of-mass (cm.) 
energy. The data on production at the lowest energy 
are consistent with the particles in the final state being all 
in relative 5*- waves, with the only feature evident in the 
measured spectra coming from the strong proton-proton 
final state interaction (FSI). The lower statistics at the 
two higher energies were sufficient to extract total cross 
sections but it was hard to draw firm conclusions regard- 
ing the differential spectra which, on general grounds, are 
expected to be much richer than at — 18.5 MeV. 

In the pioneering work of the DISTO collaboration, 
strong evidence was presented for the importance of 
higher partial waves in the pp — >■ pp(j) reaction at = 
83 MeV [13], but no attempt was made to make a con- 
sistent partial wave decomposition. The reason was in 
part due to the necessity to study in detail the struc- 
ture of the non-0 K^K~ background. We have since 
then shown [l^ that, due to the K~p final state interac- 
tion, the differential spectra for the non-0 contribution 
to the pp — ppK~^K~ reaction are strongly distorted. 
As a consequence, one needs full descriptions of both 
and non-0 components in order to extract credible par- 
tial wave parameters. When these are implemented in 
our current data set it is found that only a small amount 
of the total pp — >■ pp(j) cross section at = 76 MeV 
corresponds to pure S'-wave final states and this explains 
the non-observation of the S'-wave pp FSI enhancement 
in the results. However, the energy dependence of the to- 
tal production cross section then seems to require some 
enhancement in the (pp system at low invariant masses. 

The distortion of the K~p and K~pp invariant mass 
spectra observed in the pp — >■ ppK^K~ reaction were 
very well parametrized by assuming factorized pair- wise 
final state interactions in the pp, K~p, and K~^K~ sys- 
tems with constant effective scattering lengths [IB] ■ Fur- 
thermore, their inclusion led to a good description of the 
energy dependence of the total cross section, including 
the low energy COSY-11 data [18-20]. The distortions 
were in fact first identified in data taken below the 
threshold where selection of the non-0 contribution is 
automatic [2^. These distortions were well described by 
the FSI parameters used at higher energies but, taken 
in isolation, the error bars on the parameters extracted 



from these low energy data were very large [ll']. 

We here present much more precise differential data 
for the pp — > ppK^K^ reaction at a beam energy of 
Tp = 2.83 GeV (e^ = 76 MeV) obtained using the COSY- 
ANKE spectrometer. These will challenge the theoretical 
models that can describe well the energy dependence of 
the pp — >■ pp4> total cross sections. Currently few of the 
experimental spectra are calculated in any of the phe- 
nomenological approaches. 

The paper is organized as follows. We first describe the 
experimental setup and data analysis in Sec. |TT1 The de- 
tailed phenomenological parametrizations developed for 
and non-0 K~ production needed to make the ac- 
ceptance corrections are described here. The resulting 
differential distributions for production are presented 
in Sec IIIIl with the integrated cross sections for all the 
pp ppK^K^ data being given in Sec. IIVI The distor- 
tions in the non-0 differential cross sections arising from 
the various final state interactions are discussed in SecfVl 
followed by our conclusions in Sec. I VII 



II. EXPERIMENT AND DATA ANALYSIS 



The experiment was performed at the Cooler Syn- 
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chrotron (COSY) of the Forschungszentrum Jiilich 
using the ANKE magnetic spectrometer [2^ [^J] that is 
located at an internal target station of the storage ring. 
ANKE contains three dipole magnets; Dl and D3 di- 
vert the circulating beam onto the target and back into 
the COSY ring, respectively, while D2 is the analyzing 
magnet. There are detection systems placed to the right 
and left of the beam that register positively and nega- 
tively charged ejectiles, respectively, while fast positive 
particles are measured in the forward detector. Both 
the positive and negative side detectors consist of sets 
of start and stop scintillation counters for time-of-flight 
(TOF) measurements and two Multiwire Proportional 
Chambers (MWPCs) utilized for particle momentum re- 
construction. Two layers of scintillation hodoscopes and 
three MWPCs are incorporated in the forward detector 
which, in addition to studying a fast proton from the 
pp — > ppK^K~ reaction, is also used to measure sub- 
sidiary reactions that are needed to determine the lumi- 
nosity [111 . A high density hydrogen cluster-jet 
target l26l| was employed in the experiment. 

Particle identification relied on momentum determina- 
tion and time-of-flight measurements [H, 0, The 
time differences were calibrated by using, respectively, 
7r"*'7r~ and n^p pairs for the negative and forward STOP 
counters with respect to the positive STOP counters, de- 
scribed in detail in Ref. |24|. The kaon detection effi- 
ciency depends on the particle's momentum and varied 
between 92% and 98%, whereas that for the forward- 
going protons was about 96%. The uncertainties in the 
efficiency estimates were about 3%. 

The pp — ppK^K~ reaction was identified through 
a triple-coincidence involving the detection of a K^K~ 
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pair and a forward-going proton, with the additional re- 
quirement that the missing mass of the K~^K~p system 
be consistent with that of the proton. In the off-hne 
analysis, positive kaons were selected by a dedicated 
detection system using the TOF information between 
the START and the STOP counters The R- and 

forward-going protons were then selected simultaneously 
using the time-of-flight differences, as described in detail 
in Ref. M. 
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FIG. 1: The K K~p missing-mass distribution in the pp 
pK^ K~ X reaction at Tp — 2.83 GeV. The hatched histogram 
shows the cuts used for the selection of the non-detected pro- 
ton. The solid line, which is a second-order polynomial curve, 
was used in the analysis to estimate the background contri- 
bution under the proton peak. 



is significant and this leads to large fluctuations. Phe- 
nomenological parametrizations that describe well the 
experimental data in both the </) and non-0 regions must 
therefore be relied upon in order to perform the necessary 
acceptance corrections. 
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FIG. 2: (Color online) The raw K'^ K~ invariant mass dis- 
tribution, IMj^-+^- (points), is compared with the distribu- 
tion of events obtained in a Monte Carlo simulation (curve). 
The error bars indicate only the statistical uncertainties. The 
blue curve shows the non-(/> contributions within the fitted 
parametrization, the red curve the four-body phase-space sim- 
ulation of ppK'^ K~ , and the dotted histogram the (f) contribu- 
tions. The solid line is the incoherent sum of the (j> and non-</) 
contributions. The vertical lines indicate the cuts used for the 
separation of the cj>-nc\i and non-(j> regions. The fiuctuations 
reflect the Monte Carlo sampling effects. 



All time-of-flight selections, as well as the identification 
of the , were performed within ±3f7 ranges. A similar 
cut was also made in the missing-mass distribution of 
the detected K^K^p shown in Fig. [TJ The fraction of 
misidentified events inside the ±3ct (ct = 4.7 MeV/c^) 
cut window around the proton mass was estimated to be 
about 11.5%, which was subtracted in the analysis using 
weighted data from the side bands, as parametrized by 
the solid line. Any ambiguity in this procedure is less 
than 3% and is considered as one source of systematic 
uncertainties in the analysis. 

Having identified good pp -> ppK^K~ events, these 
were binned in terms of the K^K~ invariant mass, 
IM j^+j^-, and the corresponding results are shown in 
Fig. [2] A clear peak is observed above a slowly varying 
background. The experimental data were then divided 
into two samples, a (/)-rich region where 1.01 GeV/c^ < 
IMji+x- < 1-03 GeV/c^ and a non-(/) (the rest) region. 
The mo del- independent acceptance estimate method 
used in our earlier work [ill, 13; cannot be applied in 
the present analysis since, at this higher excess energy, 
the number of zero elements in the acceptance matrix 



We start the analysis with the kaon pair production 
away from the cj) region since it is crucial to master this 
contribution to understand the background under the (jj 
peak. The ansatz in our previous work on non-0 produc- 
tion [l^ was taken as the basis of the simulation. Here it 
was assumed that the overall enhancement factor was the 
product of enhancements in the pp and two K^p systems: 

F = Fpp{qpp) X FKpiqKpJ x FxpiqKi.^) (1) 

where qpp, qxpi and qxp-i are the magnitudes of the rela- 
tive momenta in the pp and the two K~p systems, respec- 
tively. Note that it is believed that the K^p interaction 
might be weakly repulsive and, if so, its effects would be 
interpreted as extra attraction in the K~p system. 

Using an effective K~p scattering length of a^-p — 
{0+1.5i) fm [16] , together with an additional weight of 1+ 
2.0 cos^ 9 on the polar angle of the K^K^ system in the 
overall cm. system, the invariant mass distributions can 
be described very well, except for the very low K'^K~ 
invariant masses, IM-x+k- < 995 MeV/c^. In this region 
there are small residual effects associated with KK final 
state interactions [1^, to which we shall return later. 
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Seven degrees of freedom are required to parametrize 
the unpolarized ppK~^K~ final state and these were cho- 
sen to be four angles, the K^K" and K^K^p invariant 
masses, and the relative momentum of the protons in the 
pp system. Distributions in these seven variables were 
generated inside the ANKE acceptance and compared 
with the experimental data for non-0 data in Fig. |31 
It is evident that the description of the non-resonant 
pp — ppK^K~ reaction at Tp — 2.83 GeV is very sat- 
isfactory and the same is true for the K'^ K~ invariant 
mass distribution of Fig. [5] when the (j) contribution is 
added incoherently. 

Turning now to (j) production in proton-proton colli- 
sions, the only amplitude that survives at threshold cor- 
responds to the ^Pi — >■ transition. We here denote 
the final state by ^^~^^Lj£, where S, L, and J represent 
the total spin, orbital angular momentum and total an- 
gular momentum of the pp system, respectively, and £ 
the orbital angular momentum of the relative to the pp 
system. Our previous analysis indicates that the differ- 
ential cross section at an excess energy = 18.5 MeV is 
dominantly S'-wave, with a clear effect coming from the 
pp final state interaction (l5| . 

In contrast, significant contributions from higher par- 
tial waves were suggested by the DISTO data at — 
83 MeV , where the differential cross sections as func- 
tions of the proton momentum in the pp reference frame 
and the momentum of meson in the cm. system were 
interpreted as reflecting the importance of Ps and Sp 
final waves, respectively. The anisotropy in the helicity 
distribution shows the necessity also for a Pp contribu- 
tion. 

There are several possible transitions that could lead 
to a ppcj) final state and we keep only typical ones in 
our model description, where the spin-averaged squared 
transition matrix element is written as: 



|M|2 = Ass {k X + Aps p 
+Asp hiq- kf - q 



App{q-pf 



(2) 



The momenta of the proton beam and (j) meson in the 
overall cm. system are denoted by K and q, respectively, 
k represents the momentum of decay kaons in the ref- 
erence frame, and p is the relative momentum in the final 
pp system. 

Apart from the explicit momentum factors, we assume 
that the coefficients A^i in Eq. ^ are constant except 
that, at low invariant masses, the final pp system in the 
^Sq state is subject to a very strong final state inter- 
action. The Ass and Asp contributions in Eq. ^ were 
therefore multiplied by an enhancement factor which was 
calculated using the Jost function. 



Fpp(.qpp) — \ J{.Q-P 



9i 



(3) 



pp 



where we take a = 0.1 fm and /3 = 0.5 fm ^ [l6| . 
The Coulomb interaction was neglected and, crucially. 
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FIG. 3: Differential distributions of experimental (points) 
and simulated (curves) yields for kaon pair production in the 



reaction at Tr, = 2.83 GeV for the non- 



re- 

„2\ 



pp — ^ ppK'^K 

gions (IM^+^f- < 1.01 GeV/c^ or IM^.+ _k- > 1-03 GeV/c"). 
Vertical bars represent the statistical uncertainties and hori- 
zontal ones the bin widths. The individual panels are (a) the 
cosine of the polar angle of the in the K^K~ reference 
frame, (b) the polar angle of the kaon pairs in the overall 
cm. frame, (c) the polar angle of the emitted proton in the 
pp reference frame relative to the beam direction, (d) the po- 
lar angle of the proton in the pp reference frame relative to 
the direction of the kaon pair, (e) the proton momentum in 
the pp reference frame, and (f) the K^K^p invariant mass 
distribution. 



no attempt was made to include a final state interaction 
in the (l)p system. 

A Monte Carlo phase-space simulation was written, 
based on GEANT4 [1^, which took into account the de- 
tector efficiency, resolution and kaon decay probability. 
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TABLE I: Values of the model parameters of Eq. ((2| deduced 
by comparing the simulations with data in the (j) region. The 
momenta are measured in GeV/c. All the parameters are 
normalized to Ass = 1 and the corresponding uncertainty of 
±0.25 is not included. 



Parameter 


Fit value 
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1.0 


Asp 


9.9 ±1.8 
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143 ±4 
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293 ± 21 
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cross sections for both (p and non-(/) kaon pair produc- 
tion. The luminosity needed for this analysis was deter- 
mined on the basis of the pp elastic scattering data that 
were taken in parallel, using the forward detector [l6j . 
The associated systematic uncertainty is estimated to be 
7%. Systematic uncertainties also arise from the back- 
ground subtraction, tracking efficiency, and the model- 
dependent acceptance corrections. The latter were es- 
timated from the differences between the distributions 
corrected by the parametrization and those corrected by 
phase space. As the observed distributions deviate signif- 
icantly from phase space, such estimates provide upper 
limits on these uncertainties. 



III. DIFFERENTIAL CROSS SECTIONS FOR </> 
PRODUCTION 



The angular distributions for the pp — > ppcj) reaction 
measured in this experiment and that of DISTO [l^l are 
shown in Fig. [SJ These distributions must be symmetric 
about cos 9 — Q and the data can be parametrized in the 
form: 



FIG. 4: The differential distributions of yields in the (j) re- 
gion (1.01 GeV/c^ < IMji+K- < 1-03 GeV/c^) for the 
pp — !■ ppK^K~ reaction at Tp — 2.83 GeV, where the points 
are experimental and the curves represent simulations. The 
notations for the six panels are the same as those in Fig. [3] 



— = a 1±6P2 cos( 
alt 



(4) 



The numerical values of the coefficients obtained from 
fitting the data are reported in Table [TTl 



Contributions from the phenomenological parametriza- 
tions for and non-c/) production were then included as 
weights. The (f> meson was taken to have a Breit-Wigner 
form with a width of F = 4.26 MeV/c^ [sSl, convoluted 
with a resolution width of cr ~ 1 MeV/c^. The values of 
the coefficients A^g in Eq. ([2]) were determined by min- 
imizing in the difference between the simulated and 
experimental spectra, and the results are shown in Ta- 
ble HI The resulting descriptions of the experimental data 
in Figs, m [31 and m are very good and certainly adequate 
for carrying out the acceptance corrections. 

The fits of the parametrizations to the experimental 
data allow the extraction of the differential and total 



TABLE II: Values of the coefficients of Eq. Q for the K+ 
decay angle with respect to the beam direction, the cm. pro- 
duction angle, and the helicity angle, deduced by fitting the 
data of ANKE and DISTO The DISTO data have been 
scaled by 0.7 in order to allow a direct comparison of the two 
sets of results. 





ANKE 


DISTO (scaled by 0.7) 


a [nb/sr] 


b 


a [nb/sr] 


b 


cos 

cos et.m. 
COS ^Ppp 


10.96 ± 0.23 
10.71 ± 0.21 
11.06 ±0.22 


-0.11 ±0.04 
0.21 ±0.04 
0.07 ± 0.04 


10.78 ±0.50 
10.81 ±0.45 
10.64 ± 0.64 


-0.27 ±0.08 
0.07 ±0.07 
0.30 ±0.14 
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FIG. 5: (Color online) Angular distributions of the pp — > 
PP4> reaction obtained in this experiment (solid circles), 
where the systematic uncertainties are shown by the hatched 
histograms, compared with the scaled DISTO data (open 
squares) The dotted curves represent isotropic distri- 

butions whereas the solid ones show fits to the ANKE results, 
(a) The distribution with respect to the cosine of the K'^ 
polar angle in the (j) rest frame (Decay angle). The dashed 
curve demonstrates a sin'^ 9^ behavior, (b) The distribution 
in the (p polar angle in the overall cm. system. The blue, 
red and green curves are typical theoretical predictions from 
Refs. [iJl and [l^, respectively, (c) The distribution in 
the proton polar angle in the pp reference frame relative to 
the incident proton direction (Jackson angle), (d) The dis- 
tribution of the proton polar angle in the pp reference frame 
relative to the direction (Helicity angle). 



In the near-threshold region where the Ss final state 
dominates the meson spin must lie along the beam 
direction. The polar angular distribution of the decay 
kaons in the (j) meson rest frame should then display a 
sin^ 6^ distribution, where 6^ is the angle of a daughter 
kaon from the (j) decay in the 4> rest frame. The data at 
£0 = 18.5 MeV [15^ are consistent with such a dependence 
and deviations from this behavior are a sign of higher 
partial waves. 

Quite generally, the differential cross section is of the 
form: 

^ cx [(1 - poo) sin^ 6^ + 2poo cos^ e'^] , (5) 

where poo — (1 + 6a:)/3 is a spin density matrix ele- 
ment. From the value of given in Table HI] it is seen 
that the ANKE results correspond to poo = 0-30 ± 0.01, 
which is close to the unpolarized value of i. This is to 
be compared with the value of poo = 0.23 ±0.04 reported 



by DISTO at the marginally higher = 83 MeV ^ 
where, in both cases, only statistical errors are quoted. 
These are model-independent proofs that higher partial 
waves are important at even relatively modest excess en- 
ergies. A similar conclusion is reached in a study of the 
available pn -> data [3l| . 

The angular distribution of the ip meson in the overall 
cm. frame shown in Fig. [SJd is symmetric within exper- 
imental uncertainties. The ANKE data show a stronger 
anisotropy than those of DISTO, as indicated by the 
larger b parameter in Table [III but the error bars of 
the DISTO data are significant. Theoretical calcula- 
tions [nl, [31 can describe a non-isotropy in the ex- 
perimental data reasonably well, as shown in Fig. [SJj. 
The angular distribution is expected to be isotropic when 
the mesonic current is dominant, whereas the nucleonic 
current leads to a cos^ 9 distribution. The angular distri- 
bution might therefore provide some information on the 
NN(f> coupling constant [lO, [lli • 

It could be interesting to compare our or the DISTO 
results of Table HI] with the analogous measurement at 
COSY-TOF of the pp ^ ppu reaction at an excess en- 
ergy of 92 MeV [S^ . Unfortunately, the error bars in the 
uj angular distribution, 1.0 + (0.23 ± 0.26)P2(cos 0), are 
too large to draw any useful conclusions as to whether 
the shapes are similar or not. 

The distribution in the proton polar angle measured 
in the pp reference frame relative to the beam direction 
is nearly isotropic, as shown in Fig. [5J:. This is consis- 
tent with the DISTO results. On the other hand, the 
analogous observable relative to the (f) direction shown 
in Fig. [SJi has some anisotropy. This feature, which was 
also seen in the DISTO data [17], is evidence for a con- 
tribution from a Pp final wave. 

In neither the ANKE data at = 76 MeV nor those 
of DISTO at £0 = 83 MeV is there any sign of the FSI 
enhancement in the proton-proton relative momentum 
spectrum. The lack of such an effect can be understood 
by looking at the momentum distributions of the (j) meson 
in the cm. system and relative momentum distribution of 
the final protons in the pp reference frame that are shown 
in Fig. ini The contributions of the different partial waves 
obtained by fitting Eq. ^ to the ANKE data are also 
indicated. From these it is seen that, within the given 
parametrization, the pp P-waves are completely domi- 
nant and this reduces considerably the infiuence of the 
1^0 PP FSI. 

The invariant mass distributions of the (f)p system ob- 
tained in this experiment and in the previous one at 
£0 = 18.5 MeV are presented in Fig. [71 For both en- 
ergies the data differ significantly from uniform phase- 
space predictions (dashed curve). Calculations that in- 
clude in addition the pp final state interaction (dotted 
curve) can describe the data at — 18.5 MeV, but fail 
at £0 = 76 MeV, where the higher partial waves of Eq. ([2]) 
are successful (solid curve). 
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FIG. 6: (Color online) (Upper panel) Dilferential cross section 
for the pp — !■ pp(/!> reaction as a function of the momentum of 
the (/> meson in the cm. system. (Lower panel) Differential 
cross section for the pp — >■ pp</) reaction as a function of the 
proton momentum in the pp rest frame. The systematic un- 
certainties are shown by the hatched histograms. The curves 
represent the fitted contributions of different partial waves 
within the parametrization of Eq. 



IV. THE TOTAL CROSS SECTIONS 




1.96 1.98 2.00 2.02 2.04 
(|)p invariant mass [GeV/c^] 



FIG. 7: The acceptance-corrected differential cross section 
as a function of the ^p invariant mass at excess energies (a) 
£0 = 18.5 MeV and (b) = 76 MeV. The dashed curves show 
phase-space predictions, while the dotted cures include the pp 
FSI. The solid curve represents the description of Eq. ([2]), with 
parameters being taken from Table. (T] 



sistent within statistics, though the precision of the cur- 
rent one is much higher. It should be noted that the total 
cross section for (\) production has been corrected for the 
branching ratio r^+A'-/rtot = 0.491 [33|. 



The peaking at low IM^+^f- in the raw invari- 
ant mass distribution of Fig. [5] is mainly a consequence of 
the ANKE acceptance and a smoother behavior in this 
region is seen in the acceptance-corrected data in Fig. [S] 
The contributions are there shown separately for the (p 
and non-0 contributions. Away from the low-mass re- 
gion the latter resembles quite closely that of a four-body 
ppK~^K~ phase space, which is also shown. 

The fit to the acceptance-corrected invariant mass dis- 
tribution of Fig. |5] has been used to determine separately 
the total cross sections for cj) and non-0 production mea- 
sured in the pp ppK^K^ reaction at 2.83 GeV. These 
results, together with our previous ones at this energy, 
are summarized in Table Hill The two data sets are con- 



TABLE III: Total cross section for the pp — > ppK^K~ re- 
action at Tp — 2.83 GeV separated into (p and non-0 com- 
ponents. In the case the data have been corrected for the 
— >■ K'^ K~ branching ratio. The uncertainties are, respec- 
tively, statistical and systematic. The results of previous mea- 
surements [l^, [l^ are also given. 



Channel 


cr[nb] 


a[nb] [15,^6] 


Non-0 production 
production 
Total K+K- 


91.0 ±3.0 ±11.4 
142.2 ±2.1 ±17.9 
160.8 ±3.2 ±14.4 


98.0 ± 8.0 ±15.0 
133.0 ±12.0 ±27.0 
163.3 ±10.0 ±20.0 



The total cross section for the pp — pp4> reaction is 
plotted in Fig. [H^ alongside other existing near-threshold 



8 




K"^K' invariant mass [GeV/c^] 

FIG. 8: (Color online) The acceptance-corrected pp — ^ 
ppK^K'^ dilferential cross section as a function of the K'^K~ 
invariant mass. The error bars correspond only to the statis- 
tical uncertainties; systematic uncertainties are shown by the 
hatched histograms. The blue curve shows the non-c/i con- 
tributions within the fitted parametrization, the red curve 
the four-body phase-space simulation of ppK^ K~ , and the 
dotted histogram the (f) contributions. The solid line is the 
incoherent sum of the and non-c/) contributions. 



data [m, [T3l as a function of the excess energy e^. The 
error bars shown are quadratic sums of the systematic 
and statistical uncertainties. If the coefficients A^i were 
constant, apart from the explicit momentum factors in 
Eq. ([2), then these could be used to predict the energy 
dependence of the total cross section. The resulting black 
solid curve, which by construction passes through the 
76 MeV point, underestimates severely the low energy 
data. This behaviour comes about because at 76 MeV 
the fit indicates that only a small fraction of the total 
cross section corresponds to a Ss final state and, as seen 
in Fig.[S^, the contributions from the higher partial waves 
decrease faster as threshold is approached. It therefore 
seems that there must be a strong energy variation in 
some of the A^i, which might be driven by a (pp near- 
threshold enhancement. 

The energy dependence of the total cross section is 
close to the predictions from Kaptari and Kampfer [fij, 
which include mesonic and nucleonic current contribu- 
tions. The predictions are very similar to those of three- 
body phase space modified by the effects of the pp FSI. 
This curve can fit most of the data in Fig. because, 
unlike the Ss curve, it takes the full strength at 76 MeV. 
The model of Tsushima and Nakayama also includes 
both nucleonic and mesonic current contributions but 
gives too steep an energy dependence. In neither model 




8^[MeV] 

FIG. 9: (Color online) Total cross section for the pp — > pp<j) 
reaction as a function of excess energy e^. The present result 
(solid circle) is shown together with experimental data taken 
from DISTO [l7j] (solid square) and previous ANKE measure- 
ments [l^ (open circles), a) The data are compared to the 
black solid curve derived using Eq. ([2]) with the parameters 
taken from Table. [J The individual contributions from the Ss 
(dotted curve), Sp (dashed curve), Ps (dashed-dotted curve), 
and Pp (double dotted-dashed) are normalized to their pre- 
dicted values at 76 MeV. 

b) The predictions of Tsushima and Nakayama [Tl| (ma- 
genta), scaled to pass through the 76 MeV point, underes- 
timate the low energy data. Also shown are the predictions 
of Kaptari and Kampfer ^ (green), which are very similar to 
those of three-body phase-space with the inclusion of the pp 
FSI, the results within a resonance model Xie et al. (blue), 
and a one-pion-exchange model of Sibirtsev 0] (brown), and 
this plus exotic baryons [l^ (red). 



were contributions from nucleon resonances considered 
which, if they existed, would change the energy depen- 
dence of the Ali. Also shown are the predictions of the 
resonance model of Xie et al. [13i\ . For ease of comparison, 
these have all been scaled to pass through the 76 MeV 
point. On the other hand, the one-pion-exchange calcula- 
tion , which fits the (j) production results at high energy 
(e^ > 1 GeV), fails to describe any of the near-threshold 
data. The model was subsequently extended through the 
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inclusion of baryonic resonances with masses close to the 
(j)p threshold ^J]. This achieves a better description at 
lower energies, as shown in Fig. [91d. It is clear from this 
discussion that the behavior of the total cross sections is 
insufficient by itself to distinguish between different the- 
oretical models; such calculations must be tested against 
various differential spectra of the types presented here. 

The current value of the pp — > pp(j) total cross section at 
76 MeV given in Table IIIII is only a little higher than our 
previous result [l^. The conclusion drawn there, that 
the ratio of this to the cross section for ^production is 
about a factor of six above the OZI limit 2], is therefore 
still valid. 

Values of the non-0 contribution to the pp — !■ ppK~^K~ 
total cross section were reported in our earlier work f T6j 
and any change in the 76 MeV point is well within the 
total error bars. It was shown there that the energy de- 
pendence of this cross section could only be understood 
fully if all the final state interactions in the pp, K^p, and 
K^K^ subsystems were included in the estimates. 

V. NON-(^ INVARIANT MASS DISTRIBUTIONS 

The strong K^p interaction can distort hugely both 
the K~p and K~pp invariant mass distributions, and 
this is taken into account through Eq. ([1]). The effects are 
most apparent if one forms the ratios of the differential 
cross sections in terms of the invariant masses: 

_ da/dMK~p _ da/dMji-pp 

" da/dMK^p ' ^^^^ " da/dMK+pp ' ^ ' 

The corresponding experimental data and simulations 
are to be found in Figs. [TU] and [TTJ If the K^p and 
K~p final state interactions were identical, then the ra- 
tios Rkp and Rxpp would be constant and equal to unity. 
However, both Rkp and Rkpp display very large prefer- 
ences for lower invariant masses, which probably reflect 
an attraction between the K~ and one or both of the pro- 
tons. Similar effects have been observed at lower excess 
energies [Hill HI. 

The general features of these results are well repro- 
duced by the simple factorized ansatz of Eq. ([ij. It is 
nevertheless surprising that the distortions produced by 
the constant effective scattering length, a = (0 + 1.5i) fm, 
used at ekk = 51 MeV [3 stfll describe the data so well 
at an excess energy with respect to the K^K~ threshold 
as high as 108 MeV, though some deviations are apparent 
for Kp invariant masses above about 1.5 GeV/c^. 

The distortions of Figs. [TU] and [TT] clearly indicate that 
the direct production of the scalar resonance oq or /o 
cannot be the dominant driving mechanism in the pp — >■ 
ppK~^K~ reaction. On the other hand, the strength 
of the K~'p interaction suggests that kaon pair produc- 
tion might be related to that of the A(1405) through 
pp pK+ {A{1405) K'p) [l[. This idea was put on 
a quantitative footing by assuming that the An405) was 
formed through the decay N* if+A(1405) 




1.45 1.50 

Kp invariant mass [GeV/c1 



FIG. 10: Differential cross sections for the pp — >■ ppK K~ 
reaction in the non-(j) region with respect to the invariant 
masses of K^p (upper panel) and K~p (middle panel), and 
their ratio Rkp (lower panel). The dashed histograms repre- 
sent the four-body phase-space simulations, whereas the solid 
ones represent the theoretical calculations taking into account 
pp and K~p final state interactions through Eq. ([T]). 

The simple ansatz of Eq. ([1]) underestimates the cross 
section for low K'^K^ masses, i.e., in the interval be- 
tween the K+K- and K°K° thresholds at 987.4 and 
995.3 MeV/c^, respectively. Similar effects were ob- 
served in pp ppK+K- by DISTO 17] and by ANKE 
in pn — )• dK~^K~ [s^l- Although these enhancements 
must be due to KK final state interactions, including 
K^K~ ^ K'^K'^ charge exchange scattering, they could 
be connected with some small production of the ao//o 
scalar resonances. However, in reality, the data are only 
sensitive to the KK scattering lengths. 

A combined analysis of ANKE data at three ener- 
gies [2^ suggests that, independent of the exact values of 
the scattering lengths, the KK enhancement is mainly in 



10 




2.40 



2.45 



Kpp invariant mass [GeV/c ] 



FIG. 11: Differential cross sections for the pp — > ppK'^K~ re- 
action in the non-0 region with respect to the invariant masses 
of K^pp (upper panel) and K~pp (middle panel), and their 
ratio Rkpp (lower panel). The dashed histograms represent 
the four-body phase-space simulations, whereas the solid ones 
represent the theoretical calculations taking into account pp 
and K~p final state interactions through Eq. ([1]). 



the isospin-zero channel. The model for the enhancement 
factor fitted there has been introduced into the simula- 
tion to describe better the data shown in Fig. [8] for in- 
variant masses IMkk < 995 MeV/c^. Its effects can be 
seen more clearly in the plot of the ratio of the K^K^ 
invariant- mass data to the simulation based on Eq. p]). 
where no KK FSI was included. This, together with 
the results of previous measurements |16j , are shown in 
Fig. [121 The two data sets are in agreement and are con- 
sistent with the existence of some coupled-channel effect 
at the K^K^ threshold but much better data would be 
required to prove this unambiguously. 
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FIG. 12: Ratio of the measured K~ invariant mass in 
the pp — >■ ppK^K~ reaction to estimates based on Eq. ([1}. 
In addition to the current data (solid circles), the weighted 
averages of previous measurements (open squares) 16] are 
also presented. The solid curve represents the best fit in the 
model of Ref . [28|] , which includes charge-exchange and elastic 
K^K~ FSI. The best fits neglecting charge exchange and 
including only this effect are shown by the dashed and the 
dot-dashed curve, respectively. 



VI. DISCUSSION AND CONCLUSIONS 

New measurements of the differential and total cross 
sections for the production of kaon pairs in proton- 
proton collisions have been presented at a beam energy 
of 2.83 GeV. The reaction was identified through a triple 
coincidence of a K^K~ pair and a forward-going pro- 
ton detected in the COSY-ANKE magnetic spectrome- 
ter, with an additional cut being placed on the missing- 
mass spectrum. 

By careful modeling, it was possible to describe all the 
experimental spectra in regions of the K~^K~ invariant 
mass where the </> meson sits, as well as at smaller and 
larger masses. This allowed acceptance corrections to be 
made in order to extract cross sections where the (j) and 
non-(/) contributions were reliably separated. 

The main feature of the non-cj) data is the very strong 
distortion of both the K^p and K~pp spectra by the 
K^p final state interaction. This may be a reflection of 
the excitation of the A(1405) in the production process 
and, in the Xie and Wilkin approach 33], the produc- 
tion of non-0 kaon pairs proceeds mainly through the 
excitation of i^T """-hyperon pairs. It is remarkable to note 
that these distortions are described quantitatively by the 
factorized approximation of Eq. ([T]) with the same con- 
stant scattering length that was used for the lower en- 
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ergy data [16|. On the other hand, because these data 
correspond to events where the K~^K~ emerge in the fi- 
nal state, they cannot contribute directly to the ongoing 
debate regarding the possibility of deeply bound K~pp 
states [35], except to emphasize that the K^pp interac- 
tion is still strong even above threshold. 

There is evidence for some K^K'' final state inter- 
action that changes in nature at the K^K^ threshold 
but the contribution of this region to the integrated 
pp — >■ ppK^K~ cross section is very small and it is hard 
to find any indication of the excitation of the ao//o scalar 
resonances in the reaction. As already pointed out in 
our earlier work [l^ , the energy dependence of the total 
cross section near threshold can be understood simply in 
terms of the effects of the pp, K~p, and K~^K~ FSI. To 
establish a better understanding of possible structure at 
the K^K'^ threshold, high statistics are required in this 
region and this might be achieved in the data collected 
below (j) production threshold (36| . 

Having a good description of the background, it was 
possible to derive detailed invariant mass and angular 
distributions for the pp — )■ ppcj) reaction. Although 
the DISTO collaboration [13] showed the significance of 
higher partial waves at the marginally higher excitation 
energy of — 83 MeV, they did this mainly on the 
basis of relative momentum spectra. Their conclusion 
is confirmed unambiguously by the angular distributions 
presented here. For example, at = 18.5 MeV the 0- 
meson is completely aligned, as it has to be for a Ss final 
state [15;]. In contrast, in the present data the emerging 
(j) is almost unpolarized and this clearly signals the pres- 
ence of higher partial waves. This is consistent with the 
evidence from the momentum distributions, which also 
show the dominance of P waves in the final pp system. 
This explains why the ISn pp FSI, which is so impor- 
tant at £0 = 18.5 MeV [l|, is not observed at 76 MeV. 
Furthermore, in contrast to the DISTO result flf\, clear 
anisotropy was observed in the cm. angular distribu- 
tion and this can be ascribed to the contribution from 
p wave. This angular distribution might provide infor- 
mation on nucleonic current contributions and the NN(j) 
coupling constant [l3, Uli • 

Even if one considers only a few partial waves, there 
are simply too many parameters to perform useful fits 
and only typical Ss, Sp, Ps, and Pp contributions were 
considered in Eq. ^ . The fitted data show that the con- 
tribution of the final Ss wave to the cross section repre- 
sents only a small amount of the total at = 76 MeV. 
As a consequence, the extracted parameters predict a 
total cross section that grossly underestimates the mea- 
surements at lower energies. 

The simplest way out of the total cross section dilemma 
would be to assume that a c/fp threshold enhancement 
leads to a significant energy dependence of some of the 
A^i coefficients. In this context it is interesting to 
note that the large contribution of the Pp wave to the 



pp — )■ pprj cross section at an excess energy of 72 MeV was 
ascribed to a strong T]p FSI driven by the iV*(1535) iso- 
bar [13]. Against the (j)p enhancement hypothesis is the 
fact that the large excess of events in the (pp invariant 
mass distribution shown in Fig. [7}d at low masses can be 
explained in the partial wave fitting of Eq. ([2]), without 
including any (j)p enhancement. We have not, however, 
shown that the fitting of the data is unambiguous and 
there could be other truncated partial wave forms that 
might be equally successful. Furthermore, from the start 
we have not included any final state interaction between 
(j) and protons in the parametrization. There could there- 
fore be a possible trade-off between some of the partial 
wave parameters and an FSI in the (pp system. Neverthe- 
less, the phenomenological parametrization is sufficient 
for acceptance correction and it describes well most of 
the differential distributions. 

In the parametrization of Eq. ([2|) , the coefficients A^g 
were taken to be constant and no resonance effects were 
included. Recent theoretical studies have suggested that 
bound states or resonances might be formed in the near- 
threshold (pp system 38, 39] and, if so, they would cer- 
tainly influence the behaviour of some of the A^i. In 
this context, it is interesting to note that a bump was 
observed in the near-threshold^^ meson photoproduc- 
tion from hydrogen by LEPS 4Cr] and in the prelimi- 
nary results of CLAS (41|]. Furthermore, it seems that s- 
wave production of the (p in the pd — >■ ^HeX reaction is 
anomalously large compared to the uj and 77' mesons . 
Such effects might even be part of the explanation for the 
violation of the OZI rule ^ in the ratio of to w produc- 
tion. Alternatively, it is possible that other strangeness 
production channels could influence the energy depen- 
dence of the pp -> ppK^K^ reaction psI - IZsj ]. 

Although some theoretical models have been able to 
describe a posteori the published total cross sections for (p 
production, calculations of differential distributions with 
which to compare our experimental data are rather lim- 
ited. It is only when a model is tested against a range 
of differential distributions, as presented here, that some 
credence can be given to the model. Total cross sections 
are insufficient and more theoretical work is therefore re- 
quired. 
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